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Pluripotent stem cells evade replicative senescence, whereas other primary cells lose their proliferation and differentiation 
potential after a limited number of cell divisions, and this is accompanied by specific senescence-associated DNA meth- 
ylation (SA-DNAm) changes. Here, we investigate SA-DNAm changes in mesenchymal stromal cells (MSC) upon long- 
term culture, irradiation-induced senescence, immortalization, and reprogramming into induced pluripotent stem cells 
(iPSC) using high-density HumanMethyIation450 BeadChips. SA-DNAm changes are highly reproducible and they are 
enriched in intergenic and nonpromoter regions of developmental genes. Furthermore, SA-hypomethylation in partic- 
ular appears to be associated with H3K9me3, H3K27me3, and Polycomb-group 2 target genes. We demonstrate that 
ionizing irradiation, although associated with a senescence phenotype, does not affect SA-DNAm. Furthermore, over- 
expression of the catalytic subunit of the human telomerase (TERT) or conditional immortalization with a doxycycline- 
inducible system (TERT and SV40-TAg) result in telomere extension, but do not prevent SA-DNAm. In contrast, we 
demonstrate that reprogramming into iPSC prevents almost the entire set of SA-DNAm changes. Our results indicate that 
long-term culture is associated with an epigenetically controlled process that stalls cells in a particular functional state, 
whereas irradiation-induced senescence and immortalization are not causally related to this process. Absence of SA- 
DNAm in pluripotent cells may play a central role for their escape from cellular senescence. 



[Supplemental material is available for this article.] 

Pluripotent cells have unlimited self-renewal potential. Under 
appropriate culture conditions, embryonic stem cells (ESC) as well 
as induced pluripotent stem cells (iPSC) can virtually be passaged 
infinitely without any signs of replicative senescence (Zeng 2007). 
In contrast, somatic cells enter a senescent state after a limited 
number of cell divisions (Hayflick and Moorhead 1961). The same 
holds true for the differentiated progeny of iPSC and ESC, indicat- 
ing that these cells are not immortalized per se (Lapasset et al. 2011). 

Senescent cells remain metabolically active, but their pro- 
liferation potential is unequivocally lost. They acquire typical 
morphological changes and functional deteriorations such as loss 
of differentiation potential. "Replicative senescence" is associated 
with telomere attrition during culture expansion, whereas "pre- 
mature senescence" is a consequence of other stresses such as 
oxidative stress, DNA damaging agents, irradiation, or oncogene 
activation — both of these aspects may be superimposed during in 
vitro culture. Senescence does not occur simultaneously during 
culture expansion and some cells may already have entered a se- 
nescent state, whereas others are still able to proliferate (Wagner 
2010). All of the above-mentioned mechanisms ultimately result 
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in DNA damage that triggers activation of common pathways such 
as TP53/CDKN1A and CDKN2A/RB1 (Seluanov et al. 2001; Shay 
and Roninson 2004; Bazarov et al. 2012; Sperka et al. 2012). Plu- 
ripotent cells express telomerase and maintain telomere integrity 
over time, which may prevent replicative senescence (Marion et al. 
2009b). In addition, it has been speculated that long-term culture- 
associated changes are related to a developmental process that 
might be functionally relevant for suppressing tumorigenesis in 
vivo (Campisi and d'Adda di 2007; Kang et al. 201 1). This notion is 
supported by specific senescence-associated epigenetic modifica- 
tions, pointing to a rather coordinated process. 

DNA methylation (DNAm) predominantly occurs at cyto- 
sines in the context of CpG dinucleotides (Lister et al. 2009). It 
may influence gene expression by direct interference with tran- 
scription factors or with methyl-CpG-binding proteins that mod- 
ify histones, and thereby inactivate promoter regions. Our group 
has previously demonstrated that long-term culture of mesen- 
chymal stromal cells (MSC) and fibroblasts coincides with differ- 
ential methylation at specific CpG sites (Bork et al. 2010; Koch 
et al. 2011; Schellenberg et al. 2011). These senescence-associated 
(SA) modifications are highly reproducible and continuously ac- 
quired in the course of culture expansion — they can even be used 
as biomarkers to account for the number of passages or the time 
of in vitro culture (Koch et al. 2012). 

In this study we have further analyzed SA-DNAm changes. 
MSC were subjected to (1) long-term culture, (2) irradiation-induced 
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senescence, (3) immortalization, and (4) reprogramming into 
induced pluripotent stem cells, which we refer to as iP-MSC. 
Subsequently, DNAm was analyzed using a novel Infinium 
HumanMethylation450 platform that can assay more than 480,000 
CpG sites at single-base resolution (covering 99% of RefSeq genes 
and 96% of CpG islands) (Bibikova et al. 2011). We demonstrate 
that SA-DNAm changes are significantly enriched in developmental 
genes such as homeobox genes. Particularly, SA-hypomethylation 
occurs in intergenic regions or nonexpressed genes. Neither se- 
nescence induction by gamma irradia- 
tion, nor telomere extension influence 
SA-DNAm. In contrast, the majority of « 
SA-DNAm is prevented by reprogram- 
ming into iP-MSC. Therefore, SA-DNAm 
changes do not resemble stochastic mod- 
ifications that accumulate over time — 
they rather reflect a tightly regulated epi- 
genetic program that is prevented by 
reprogramming into a pluripotent state. ^ 



of them were restricted to circumscribed genomic locations. CpG 
sites were classified according to gene regions as described before 
(Sandoval et al. 2011; Supplemental Fig. 2). Unexpectedly, SA- 
CpGs were hardly located within 200 bp upstream of the tran- 
scription start site (TSS200; P< 10~ 18 ; hypergeometric distribution) 
and in the first exon (P < 10~ 5 ) — those regions where differential 
methylation has a major impact on gene expression. In contrast, 
SA-hypermethylation was particularly enriched in the gene body 
(P < 10" 7 ), 3'UTR (P < 10" 5 ), and intergenic regions (P < 10" 4 ), 



Results 

Long-term culture-associated DNA 
methylation changes 

Mesenchymal stromal cells from human 
bone marrow were expanded until they 
entered growth arrest after 80 ± 25 d and 
35.1 ± 5.7 cumulative population dou- 
blings (cPD) (Fig. 1A). At early passage, all 
cell preparations fulfilled the commonly 
used criteria for definition of MSC such as 
fibroblastoid growth, surface-marker ex- 
pression (CD14 - , CD29 + , CD31", CD34" 
CD45", CD73 + , CD90 + , and CD105 + ), 
and differentiation potential toward os- 
teogenic, adipogenic, and chondrogenic 
lineages (Dominici et al. 2006). Senescent 
MSC acquired large and flat morphology 
and expressed senescence-associated beta- 
galactosidase (SA-^-gal). The immunophe- 
notype was maintained throughout cul- 
ture expansion, but MSC of late passages 
lost their adipogenic differentiation po- 
tential and revealed reduced osteogenic 
differentiation potential (Fig. IB; Supple- 
mental Fig. 1A-D). DNAm profiles were 
compared between MSC of early passage 
(P2 or P3; n = 5) and MSC of late passage 
just before entering senescence (P7-P16 
according to inter-individual variation in 
long-term growth curves; n = 4). To focus 
on the most relevant SA-DNAm changes, 
we have only considered CpG sites with 
at least 20% differential methylation and 
a limma-adjusted P- value of <0.05. Under 
these stringent criteria 1702 CpG sites 
were SA-hypermethylated, whereas 2116 
CpG sites were SA-hypomethylated (Fig. 
1C; Supplemental Tables 1,2). 

DNAm changes were not simulta- 
neously acquired within genes, but most 
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Figure 1. (Legend on next page) 
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whereas SA-hypomethylation was highly enriched in intergenic 
regions (P < 10" 55 ) (Fig. ID). Analysis of SA-CpG sites in the con- 
text of CpG islands (CGI) revealed preferential occurrence of 
SA-hypermethylation in the 2-kb flanking regions of CGIs (shore- 
region; P < 10~ 5 ), whereas SA-hypomethylation hardly occurred 
in the context of CGIs at all (P < 1(T 264 ) (Fig. IE). SA-CpG sites did 
not exhibit clear association with subtelomeric or centromeric 
regions. We searched for differentially methylated regions (DMR) 
where 100 adjacent CpG sites on the BeadChip are regulated 
accordingly — only two genomic regions revealed such a highly 
significant enrichment of SA-DNAm changes (P < 10~ 5 ): the 
HOXB locus on chromosome 17 (SA-hypermethylated) and the 
keratin associated protein (KRTAP) locus on chromosome 21 (SA- 
hypomethylated; Supplemental Fig. 3). Overall, spatio-temporal 
DNAm changes were highly reproducible in independent samples. 

Gene expression profiles from MSC of early passage (P3-P5; 
n = 3) and late passage (P10; n = 3) were compared by microarray 
technology. Signal intensity was used as an indicator for gene ex- 
pression level and correlated with DNAm profiles. As expected, 
higher DNAm level within 200 nucleotides upstream of the tran- 
scription start site (TSS200) and first exon coincided with markedly 
reduced gene expression levels. However, the specific DNAm 
changes upon long-term culture were hardly reflected by differ- 
ential gene expression (Supplemental Fig. 4) and similar results 
were observed when we matched SA-DNAm with our previously 
published gene expression data on long-term culture of MSC 
(Wagner et al. 2008). Overall, SA-hypomethylated CpG sites re- 
vealed a mean signal intensity that was significantly below the 
average of all genes (Fig. IF). Thus, most of these genes seem to be 
neither expressed at early nor at late passages. Gene ontology 
analysis of genes with SA-hypermethylated CpG sites or SA- 
hypomethylated CpG sites revealed highly significant enrich- 
ment of genes in categories associated with cellular differentia- 
tion (Supplemental Table 3). Very similar results were observed 
when we focused on genes with SA-CpG sites within the pro- 
moter region (TSS1500, TSS200, and 5'UTR). Using gene set 
enrichment analysis we found the most significant overlap of 
either SA-hypermethylated or SA-hypomethylated CpG sites 
within genes harboring trimethylated lysine 9 of histone H3 
(H3K9me3; P = 2.6 X 10" 13 and 8.2 X 10" 13 , respectively), and 
targets of Polycomb-group family members EED (P = 2.1 X 10~ 9 
and 1.6 X 10" 7 ) and SUZ12 (P = 9.3 X 10" 8 and 2.6 X 10" 7 ) in 
ChIP array data of human ESC (Lee et al. 2006; Ben-Porath et al. 
2008). Furthermore, SA-hypomethylated CpG sites were also 
enriched in H3K27 targets of various cancer cells (Schlesinger et al. 
2007), P = 4.0 X 10" 6 ; (Kondo et al. 2008), P= 2.4 X 10" 5 ; (Acevedo 



Figure 1. Long-term culture-associated changes in MSC. (A) Long-term growth curves demonstrate 
that MSC reach a senescent state after 2-3 mo of culture expansion. (B) Cells at higher passage acquire 
morphological changes, senescence-associated beta-galactosidase staining, and loss of adipogenic 
differentiation potential (green: BODIPY-staining of lipid droplets; blue: DAPI) and decline in osteogenic 
differentiation potential (Alizarin Red staining of calcium phosphate precipitates; size bars 1 00 |xm). (C) 
Comparison of DNAm between early (P2-P3) and late (P7-P1 6) passages. CpG sites with relevant SA- 
hypermethylation (red: 1 702 CpG sites) and SA-hypomethylation (green: 21 1 6 CpG sites) are depicted 
(adjusted P-value < 0.05 and DNAm change >20%). (D) SA-DNAm changes were classified according to 
gene regions, and (f) in relation to CpG islands. (F) Signal intensity of gene expression microarrays was 
used as an indicator for gene expression level. In particular, genes with SA-hypomethylated CpG sites 
were significantly less-expressed than the average. (C) Genes with SA-DNAm changes were sub- 
sequently mapped to genes with H3K9me3, H3K27me3 (Delbarre et al. 201 0), and PRC2 targets (Wei 
et al. 2011) in human MSC. Enrichment analysis revealed that particularly genes with SA-hypometh- 
ylated CpG sites have been associated with these repressive histone marks. As an alternative, we focused 
on the subset of promoter-associated CpG sites (TSS1500, TSS200, or 5'UTR; striped bars), and this 
analysis also revealed association of genes with SA-hypomethylated CpG sites with those genes that 
have repressive histone marks (H3K9, P< 10~ 28 ; H3K27, P< 10~ 13 ; EZH2, P< 10~ 2 ). 



et al. 2008), P = 5.1 X 10" 4 . Enrichment of SA-hypermethylation 
or SA-hypomethylation was further analyzed with published ChIP 
array data for H3K9me3 and H3K27me3 in human MSC from adi- 
pose tissue (Delbarre et al. 2010) and for the Polycomb repressive 
complex 2 (PRC2) targets in human MSC from bone marrow (Wei 
et al. 2011); these marks were highly significantly enriched in 
genes with SA-hypomethylated CpG sites (Fig. 1G). Taken to- 
gether, particularly SA-hypomethylation preferentially occurred 
in regions with repressive H3K9me3 and H3K27me3 histone 
marks and in PRC2 targets — this may contribute to the very low 
expression of the corresponding genes. 

Gamma-irradiation does not affect DNAm 

Ionizing radiation (IR) of MSC induced very similar phenotypic 
changes as replicative senescence. Treatment with 15 Gy com- 
pletely abrogated proliferation and the majority of cells stayed 
metabolically active for more than one month (Fig. 2A). After 7 d 
the cells acquired the typical senescent morphology. Furthermore, 
IR impaired differentiation toward adipogenic and osteogenic line- 
age as described above (Fig. 2B; Supplemental Fig. 5), transiently up- 
regulated CDKN1A, CDKN2A, and TPS3 (Supplemental Fig. 6) and 
other investigators demonstrated activation of the corresponding 
signaling cascades on a protein level (Cmielova et al. 2012). DNAm 
profiles were compared before and 7 d after treatment with 15 Gy 
(n = 4). None of the CpG sites revealed significant changes upon ir- 
radiation (limma-adjusted P < 0.05) (Fig. 2C). Accordingly, radiation- 
induced senescence does not seem to be associated with SA-DNAm 
changes. Subsequently, we measured telomere length by mono- 
chrome multiplex quantitative PCR (Cawthon 2009). As expected, 
telomere length decreased by long-term culture (P = 0.01), whereas 
irradiation did not affect telomere length. Despite phenotypic simi- 
larity of long-term culture-associated senescence and irradiation- 
induced senescence, the causative mechanism seems to be different. 

Immortalization does not reverse SA-DNAm 

Overexpression of TERT, the catalytic subunit of human telo- 
merase, is known to extend the lifespan of primary cells, but it is 
generally insufficient for immortalization (Bodnar et al. 1998; 
Kiyono et al. 1998). We retrovirally transfected MSC with a TERT- 
IRES-GFP construct in comparison to GFP-control. Constitutive 
overexpression resulted in growth advantage of TEPT-transfected cells 
(Fig. 3A), but this was only observed in one of five experiments, and 
the cells entered a senescent state after 5 additional cPD in com- 
parison to GPP-control (Supplemental Fig. 7). Telomere length was 
moderately increased by TERT over- 

expression, whereas SA-DNAm changes 

were not affected (Fig. 3B,E). 

Alternatively, we immortalized MSC 
with a lentiviral tetracycline-inducible ex- 
pression system for TERT and SV40 large 
T antigen (SV40-TAg). In the presence of 
doxycycline (Dox) these cells could be 
passaged for more than 80 passages with- 
out any signs of senescence. However, re- 
moval of Dox from the culture medium 
resulted in immediate proliferation arrest, 
and within 10 d the cells acquired a typi- 
cal senescent phenotype with SA-(3-gal 
expression (Fig. 3C; Supplemental Fig. 8). 
qRT-PCR analysis revealed that TERT was 
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Figure 2. Irradiation-induced senescence does not result in SA-DNAm changes. (A) Representative 
growth curves of MSC after irradiation with different doses at passage 2 (0-25 Gy; black arrow). 

(B) Upon irradiation with 1 5 Gy, MSC acquire the typical senescent morphology with proliferation arrest 
and moderate SA-^-gal staining. Furthermore, adipogenic (green: BODIPY; blue: DAPI) and osteogenic 
differentiation potential (red: Alizarin Red; size bars 100 |xm) is impaired after irradiation with 15 Gy. 

(C) Scatterplot analysis of DNAm of nonirradiated versus irradiated MSC. Long-term culture-associated 
hypermethylated (red) and hypomethylated CpG sites (green) are depicted. Irradiation did not result 
in any significant DNAm changes. 



highly expressed, whereas it was hardly detectable in non- 
transfected MSC (Supplemental Fig. 6). Telomere length was 
massively extended between passages 13 and 62 (Fig. 3E). We 
used our recently published epigenetic-senescence-signature to 
test this conditionally immortalized cell line for SA-DNAm changes 
(Koch et al. 2012). DNA was harvested at passage 13 and passage 
62 (with and without Dox for 10 d) and DNAm was analyzed by 
pyrosequencing of the six relevant CpG sites. The predictions for 
passage numbers were always related to the real passage number 
(Fig. 3D; Supplemental Fig. 9). Taken together, neither IR nor TERT 
overexpression or immortalization with TERT/SV40-TAg interfered 
with the SA-DNAm changes. 

Reprogramming into iP-MSC prevents SA-DNAm 

Subsequently, we analyzed SA-DNAm of induced pluripotent MSC 
(iP-MSC). MSC at passage 3 were retrovirally infected to express 
POU5F1 (OCT4), SOX2, MYC, and KLF4 (Takahashi et al. 2007). 
After 3-4 wk colonies with typical ESC-like morphology emerged, 
expressing pluripotency markers such as POU5F1, NANOG, and 
TRA-1-80 (Fig. 4A), whereas the expression of retrovirally encoded 
reprogramming factors was silenced and PluriTest analysis of 
global gene expression profiles supported the notion that iP-MSC 
were pluripotent (Muller et al. 2011; Shao et al. 2012). Furthermore, 
embryoid bodies derived from iP-MSC comprised derivatives of all 
three embryonic germ layers. For comparison we utilized three 
widely used human ESC lines (HI, H9, and HES2). Telomere length 
differed markedly between these three ESC lines, whereas iP-MSC 



always revealed significant telomere ex- 
tension as compared with their parental 
MSC (P = 0.0018) (Fig. 3E). 

DNAm profiles were analyzed in 
iP-MSC that have been cultured for a 
similar number of passages as senescent 
MSC (P9-P11). Reprogramming of iP- 
MSC was associated with far-reaching 
epigenetic alterations and a global in- 
crease in DNAm level (Supplemental 
Fig. 10). Unsupervised hierarchical clus- 
ter analysis clearly segregated pluripo- 
tent from nonpluripotent cells and dem- 
onstrated the close relationship of iP-MSC 
and ESC. Notably, MSC of early and late 
passage were also separated (Fig. 4B). Most 
of the SA-DNAm changes were prevented 
by induced pluripotency. As expected, 
this association was more pronounced 
when we compared iP-MSC to MSC of late 
passage rather than MSC of early passage 
(Fig. 4C): 95% of the SA-hypomethylated 
CpG sites were significantly higher meth- 
ylated in iP-MSC (2027 of 2116 CpG sites; 
P < 0.05). Conversely, 46% of the SA- 
hypermethylated CpG sites were signif- 
icantly less methylated in iP-MSC in 
comparison to late passage MSC (796 of 
1702 CpG sites) (Fig. 4D). Notably, the 
same tendency was also observed when 
we compared iP-MSC with MSC of early 
passage that have so far acquired rela- 
tively few SA-DNAm changes: 87% of 
the SA-hypomethylated CpG sites re- 
vealed a significantly higher methylation in iP-MSC (1842 of 
2116 CpG sites), whereas 25% of the SA-hypermethylated CpG 
sites were significantly hypomethylated (419 of 1 702 CpG sites; 
Supplemental Fig. 11). Comparison of DNAm profiles of MSC 
and iP-MSC may be difficult because of the stark differences 
between the two cell types. Yet, it is remarkable that pluripotent 
stem cells hardly revealed SA-DNAm changes. 

We reasoned that those SA-DNAm changes particularly 
might be relevant, which appear to be prevented by pluripo- 
tency. SA-CpG sites were filtered to be differentially methylated 
in MSC of late passage in comparison to (1) MSC of early passage, 
(2) iP-MSC, and (3) ESC (adjusted P-value < 0.05; DNAm change 
>20%). A total of 404 CpG sites were overlappingly hyper- 
methylated in MSC of late passage (Fig. 4E), which were signifi- 
cantly enriched in genes with sequence-specific DNA-binding 
transcription factor activity (GO:3700; P < 10~ 12 ), particularly 
those with a homeobox domain (Pfam, P < 10~ 6 ). They were also 
enriched in genes with promoter regions, which comprise mo- 
tives for binding of POU2F1 (P = 1.1 X 10" 4 ), growth factor in- 
dependent 1 transcription repressor (GFI1; P = 1.4 X 10~ 4 ), and 
GATA4 (P = 3.5 X 10~ 4 ). Furthermore, they showed enrichment 
in targets of MIR323B (P = 2.1 X 10" 5 ) and MIR10A/MIR10B (P = 
4.6 X 10~ 4 ). In contrast, 1930 CpG sites were overlappingly 
hypomethylated at late passage (Fig. 4F), and these sites were 
significantly enriched in G-protein-coupled receptors (GO:4930; 
P = 1 .6 X 10~ 8 ) and in genes with binding sites for POU6F1, a class 
6 homeobox protein (P = 0.0025), STAT5A (P = 0.016), and GATA2 
(P = 0.017). 
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When we used our epigenetic-senescence-signature, the 
predicted passage numbers for iP-MSC and ESC were close to 
noncultured cells, although they have been culture-expanded 
over many passages (Fig. 3D). Several other investigators have 
previously compared DNAm changes upon reprogramming (Doi 
et al. 2009; Lister et al. 2009). Publicly available data sets with 
InfiniumHumanMethylation27 BeadChip were retrieved and 
tested with the Epigenetic-Sensecence-Signature (Harness et al. 
2011; Nishino et al. 2011; Nazor et al. 2012)— these pluripotent 
samples were always estimated close to noncultured cells (Sup- 
plemental Fig. 12). 

Discussion 

SA-DNAm changes are tightly coordinated and may entail func- 
tional consequences of cellular senescence. On the other hand, 
SA-DNAm does not seem to be the immediate trigger of cellular 
senescence. Irradiation-induced senescence, although phenotypi- 
cally very similar to replicative senescence, was not associated with 
SA-DNAm changes. Furthermore, SA-DNAm was not influenced 
by TERT overexpression or immortalization. We demonstrate for 
the first time that most of the SA-DNAm changes do not occur 
in iP-MSC — and this supports the notion that these epigenetic 
modifications are functionally relevant. Our results indicate that 
long-term culture is associated with a process of cellular differen- 
tiation, which continuously locks specific genomic regions in an 
inactive state. Prevention of SA-DNAm upon reprogramming into 
iPSC may be crucial for rejuvenation and escape from senescence 
(Fig. 5). 

MSC resemble heterogeneous cell preparations — the compo- 
sition of subpopulations is greatly affected by culture conditions 
and changes during culture expansion (Wagner and Ho 2007; 
Schellenberg et al. 2012). Therefore, future work should aim for 
analysis of SA-DNAm changes on single cell level. On the other 
hand, similar SA-DNAm changes were also observed in dermal fi- 
broblasts (Koch et al. 2011, 2012), indicating that the observations 
made in this study are not restricted to MSC. 

We describe how SA-CpG sites are enriched in the gene body, 
3'UTR, and intergenic regions. CpG-methylation in these regions 
is not directly related to gene expression levels (Ball et al. 2009). 
Particularly, SA-hypomethylated CpG sites are associated with 
genes showing neither expression at early, nor at late passage. In 
fact, none of the genes of the Epigenetic-Senescence-Signature 
was detected by qRT-PCR (Koch et al. 2012). This might explain 
why overall the SA-DNAm did not correlate with differential gene 
expression. Even though, SA-DNAm changes may be functionally 
relevant to permanently lock genes in an inactive state. DNAm 
directly results in tighter wrapping of DNA around the histone core 
(Lee and Lee 2012) and it may support reorganization of the nu- 
clear compartment by methyl-CpG-binding proteins (MBPs) (Sasai 



and Defossez 2009). It has been suggested that DNAm is not the 
primary determinant of gene activation in MSC but rather post- 
translational modifications of histones (Collas 2010; Sorensen 
et al. 2010). DNAm may feedback on the histone code — potentially 
by changing the binding affinity of noncoding RNAs. PRC2 is 
recruited to target genes by noncoding RNAs such as HOTAIR, 
a transcript from the HOXC locus (Rinn et al. 2007); Xist, a non- 
coding RNA involved in X-chromosome inactivation (Zhao et al. 
2008); and by the antisense-RNA 1 of CDKN2B (CDKN2B-AS1) (Yap 
et al. 2010). A recent comparison of epigenomic landscapes be- 
tween human ESC and fibroblasts revealed that the repressive 
H3K9m3 and H3K27me3 marks expand upon differentiation into 
fibroblasts (Hawkins et al. 2010). In this study, we demonstrate 
that SA-hypomethylation, in particular, is significantly enriched 
in regions with H3K9me3, H3K27me3, and targets of PRC2, and is 
in line with previous observations (Schellenberg et al. 2011). It has 
been shown that aging induces DNAm in PRC2 target genes and 
may thereby predispose to cancer (Teschendorf f et al. 2010). Many 
cancers are derived from epithelial cells, and our previous work 
indicated that similar SA-DNAm changes are also observed in other 
cell types including epithelial cells (Koch et al. 2012). This is of 
particular relevance, as SA-DNAm changes might be involved in 
malignant transformation. In fact, hypomethylation of CpG sites 
that are highly methylated in ESC has recently been shown to be a 
key determinant of tumor invasion and clinical outcome (Zhuang 
et al. 2012). 

The irreversible nature of senescent cells is generally believed 
to be due to alterations in the chromatin structure. Upon senes- 
cence, chromatin undergoes dramatic remodeling through the 
formation of facultative heterochromatin domains called senes- 
cence-associated heterochromatic foci (SAHF) (Narita et al. 2003, 
2006). These compacted foci of DNA disappear only 7 d after 
reprogramming of senescent cells, and this was even considered as 
a first sign of iPSC generation (Lapasset et al. 2011). Further re- 
search will be necessary to demonstrate whether the SA-DNAm 
reported here contributes to the formation of SAHF. 

Ionizing radiation has been shown to induce a senescent 
phenotype in various different cell types (Di Leonardo et al. 1994; 
Igarashi et al. 2007); yet it is unclear whether the underlying mo- 
lecular mechanisms are related to long-term culture-associated 
senescence (Rando and Chang 2012). We demonstrate that irra- 
diation of MSC results in permanent proliferation arrest, typical 
senescent morphology, and loss of adipogenic and osteogenic dif- 
ferentiation potential. In contrast to replicative senescence, this 
stress-induced senescence was not associated with telomere attri- 
tion and did not evoke reproducible SA-DNAm changes 7 d after 
irradiation. It is conceivable that epigenetic alterations are initially 
not present at the level of DNAm, but rather on other levels of 
the epigenome: It has been suggested that factors involved in 
transcriptional repression and gene silencing are recruited to the 



Figure 3. Immortalization does not affect SA-DNAm. (A) MSC at passage 2 were either retrovirally transfected with an TERT-\RES-GFP construct or with 
IRES-GFP (control; n = 5). Flow cytometic analysis of one sample reveals a growth advantage of 7f/?7-transduced cells at late passages. (fi)The DNAm profile 
of this 7f/?7-transfected sample was compared with the corresponding control (both passage 14). CpG sites that were hypermethylated (red) or hypo- 
methylated upon long-term culture (green) are indicated. Overall, TERT expression does not affect SA-DNAm changes. (C) MSC were immortalized by 
tetracycline-inducible expression of TERT and SV40-TAg. Upon Dox removal from the culture medium, these cells acquire a typical senescent morphology 
within 1 0 d, express SA-^-gal, and lose osteogenic differentiation potential. (D) The Epigenetic-Senescence-Signature was used to estimate cellular aging 
of all samples (Koch et al. 201 2). DNAm was determined at six SA-CpG sites by pyrosequencing and used for the linear regression model to predict the 
number of passages. The predictions correlate with the real passage number for all MSC samples, including irradiated, 7f/?7-transfected, and TERT/SV40- 
TAg immortalized cells. In contrast, senescence-predictions for iP-MSC and ESC rather correspond to noncultured cells, (f) Relative telomere length was 
determined by the telomere to single-copy gene (T/S) ratio with quantitative multiplex PCR (each sample measured in triplicate). Telomere length is 
significantly reduced in MSC of later passage, whereas it is not affected by irradiation. TERT expression, TERT/SV40-TAg immortalization (particularly at 
passage 62), and generation of iP-MSC result in telomere extension. 
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Figure 4. Re prog ramming into iP-MSC prevents SA-DNAm changes. (A) Induced pluripotent MSC 
reveal typical ESC-like morphology and express the pluripotency markers POU5F1, NANOG, and 
TRA-1-80. (B) Unsupervised hierarchical clustering of global DNAm profiles clearly separates pluripo- 
tent cells (iP-MSC and ESC) and it discerns MSC of early and late passage. (C) Scatter plot of DNAm 
profiles of iP-MSC in comparison to MSC of late passage (P7-P16). In particular, CpG-sites that are 
hypomethylated upon replicative senescence (green) are hypermethylated in iP-MSC. In contrast, SA- 
hypermethylated CpG sites (red) are rather hypomethylated in iP-MSC. (D) Differential DNAm upon 
long-term culture is plotted against methylation changes upon reprogramming into iP-MSC. Overall, 
SA-DNAm changes are prevented in iP-MSC. (f,F) Venn diagrams demonstrating the overlap of SA- 
DNAm changes with differential methylation between MSC of late passage with either iP-MSC or ESC. 



Immortalization is often considered 
a counterpart of cellular senescence, but 
we demonstrate that it does not reverse 
SA-DNAm changes. As described by many 
other studies, overexpression of TERT 
may provide a proliferative advantage, 
but does not assure immortalization of 
MSC (Bodnar et al. 1998; Kiyono et al. 
1998). In contrast, combination of TERT 
with SV40-TAg facilitates reliable immor- 
talization (O'Hare et al. 2001). In this 
study, we describe for the first time a 
conditionally immortalized MSC line with 
TERT and SV40-TAg using a tetracycline- 
inducible system. Despite immortaliza- 
tion, the in vitro differentiation poten- 
tial decayed upon long-term culture — 
adipogenic differentiation could not be 
induced at passage 13 and osteogenic 
differentiation decreased at higher pas- 
sages. TERT overexpression, as well as 
TERT/SV40-TAg immortalization, resulted 
in extended telomeres, but they did not 
affect SA-DNAm. Therefore, immortali- 
zation of MSC does not address the regu- 
latory mechanisms of SA-DNAm, and 
these cells may still undergo functional 
changes throughout culture expansion. 

Telomere homeostasis is also crucial 
for pluripotent cells. Telomere shorten- 
ing is a banier for iPSC generation (Marion 
et al. 2009b), whereas transfection of 
TERT and SV40-TAg in combination with 
the classical four factors (POU5F1, SOX2, 
MYC, and KLF4) has been shown to en- 
hance the reprogramming efficiency (Park 
et al. 2008). Telomere extension upon 
reprogramming may be due to telome- 
rase expression, or it may be facilitated 
by alternative lengthening of telomeres 
(ALT), which is based on homologous re- 
combination between telomeric sequences 
(Dunham et al. 2000; Lapasset et al. 
2011). Epigenetic modifications such as 
H3K9me3, H4K20me3, and histone acet- 
ylation of telomeric chromatin have been 
shown to be involved in the regulation 
of telomere length (Benetti et al. 2007; 
Schoeftner and Blasco 2009). Furthermore, 
a decrease in DNAm, particularly at sub- 
telomeric regions, is accompanied by dra- 
matically elongated telomeres (Gonzalo 
et al. 2006). SA-DNAm was not enriched 
in subtelomeric regions, and it is yet un- 



site of DNA damage, including H3K9me3, H3K27me3, and re- 
pressive Polycomb complexes (O'Hagan et al. 2008; Chou et al. 
2010; Seiler et al. 2011). This gene silencing may ultimately acti- 
vate the same cell-intrinsic senescence program as described for 
cellular senescence — but the phenotypic and functional changes 
observed 7 d after irradiation are not causally related to DNAm 
changes at specific sites in the genome. 



clear whether or not it influences telo- 
mere length. To this end, our results demonstrate that telomere 
homeostasis is not the driving mechanism of SA-DNAm changes. 

Replicative senescence has been described as a barrier to re- 
programming of iPSC as the reprogramming efficiency is much 
higher in cells of early passage. A series of publications indicated 
that particularly up-regulation of TP53, CDKN1A, and CDKN2A 
may be responsible for this loss of reprogramming efficiency 
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Figure 5. Schematic overview of how SA-DNAm changes may influence cellular senescence. Senes- 
cence is associated with a continuous gain of SA-DNAm changes that are enriched in genes with re- 
pressive histone marks. They may be associated with permanent gene silencing, but this link is yet 
unproven and therefore labeled with a question mark. SA-DNAm changes are not affected by ionizing 
radiation or telomere length. However, these processes evoke chromosomal aberrations and may ul- 
timately activate the same downstream signal cascades of cellular senescence (including activation 
of CDKN1A, CDKN2A, and TP53). Pluripotent cells hardly reveal SA-DNAm changes; avoidance 
from these epigenetic modifications might contribute to the escape of pluripotent cells from cellular 
senescence. 



(Banito et al. 2009; Kawamura et al. 2009; Li et al. 2009; Marion 
et al. 2009a; Utikal et al. 2009). These factors are up-regulated 
upon senescence, and their ablation increased the reprogram- 
ming efficiency. It is also conceivable that the donor age re- 
sembles a barrier to reprogramming — other investigators dem- 
onstrated that it is possible to generate iPSC from centenarians 
(Lapasset et al. 2011), and it will be of interest to systematically 
analyze whether epigenetic rejuvenation is hampered in such 
samples. Overall, counteraction of senescence-associated modi- 
fications seems to play a central role for the generation of iPSC. In 
this regard, SA-DNAm might contribute to the reprogramming 
barrier at later passages. 

Conclusion 

All signs of cellular aging seem to be removed upon reprogram- 
ming into iPSC. It has even been demonstrated that continuous 
passaging of iPSC further diminishes differences in DNAm be- 
tween iPSC and ESC (Nishino et al. 2011). Yet, upon rediffer- 
entiation into fibroblasts they reveal typical signs of cellular aging, 
and finally enter a senescent state (Lapasset et al. 2011). Our 
analysis reveals that SA-DNAm changes are tightly regulated. They 
are significantly enriched in genes with H3K9me3, H3K27me3, 
and PRC2-binding sites, indicating that SA-DNAm changes are 
associated with heterochromatin formation. Overall, our results 
support the notion that cellular aging is not only due to accumu- 
lation of cellular defects, but rather resembles an organized process 
that stalls cells in a particular functional state. It remains to be 
demonstrated whether SA-DNAm changes are directly related to 
proliferation (replicative senescence) or rather some other type of 
culture stress (premature senescence). Furthermore, it is still not 
proven whether SA-DNAm resembles the initial cause of senes- 
cence or whether it is a mere marker. MSC are heterogeneous — also 
with regard to cellular aging. It is therefore conceivable that our 
iP-MSC clones were derived from a subset with little SA-DNAm, 
and that this low level is then maintained throughout culture 
expansion. Alternatively, SA-DNAm might even be reversed in the 



course of reprogramming; this would 
suggest that pluripotency entails epige- 
netic rejuvenation of cellular senescence. 

Methods 

Expansion of mesenchymal stromal 
cells 

MSC were isolated from tibia plateau of 
patients undergoing orthopaedic surgery 
(age range: 31-73 yr). All samples were 
taken after written consent using guide- 
lines approved by the Ethic Committee 
on the Use of Human Subjects at the 
University of Aachen (Permit number: 
EK128/09). MSC were culture-expanded 
in standard medium consisting of DMEM 
(PAA; 1 g/L glucose) supplemented with 
glutamine (PAA), penicillin/streptomycin 
(PAA), and 10% human platelet lysate 
(hPL) as previously described (Horn et al. 
2010). Cells grew at 37°C, in 5% CO z , 
and atmospheric oxygen conditions, and 
they were harvested by trypsinization 
upon 80% confluent growth, counted 
with a Neubauer chamber (Brand, Wertheim, Germany), and 
reseeded at a density of 5000 cells/cm 2 . cPD were calculated as 
described before (Cholewa et al. 2011). 



Immunophenotypic analysis 

Surface marker expression was analyzed on a FACS canto II (Becton 
Dickinson & Co.) upon staining with the following antibodies as 
described before (Wagner et al. 2008): CD14-allophycocyanin (APC, 
clone M5E2, BD), CD29-phycoerythrin (PE, clone MAR4, BD), 
CD31-PE (clone WM59, BD), CD34-APC (clone 8G12, BD), CD45- 
APC (clone HI30, BD), CD73-PE (clone AD2, BD), CD90-APC (clone 
5E10, BD), CD105-fluorescein isothiocyanate (FITC, clone MEM- 
226, ImmunoTools). 

In vitro differentiation 

Adipogenic and osteogenic differentiation of MSC was induced 
as described before (Pittenger et al. 1999). After 18 d, osteogenic 
differentiation was analyzed by Alizarin Red staining. Adipo- 
genic differentiated cells were stained with the green fluores- 
cent dye BODIPY (4,4-difluoro-l,2,5,7,8-pentamethyl-4-bora-3a, 
4a-diaza-s-indacene) and counterstained with DAPI (4',6-Diamidin- 
2-phenylindol; both Molecular Probes). Quantification was 
performed as described in our previous work (Koch et al. 2011). 
Chondrogenic differentiation was induced in micromass cul- 
ture for 18 d. Subsequently, the pellets were fixed with 10% 
formalin and first embedded in agarose and secondly in paraf- 
fin; 1 mm sections were stained as previously described (Koch 
etal. 2011). 



Senescence-associated beta-galactosidase 

Expression of pH-dependent senescence-associated beta-galactosi- 
dase (SA-p-gal) was analyzed using the SA-(3-gal staining kit (Cell 
Signaling Technology) or by flow cytometry with the fluorogenic 
substrate 5-dodecanoylaminorluorescein di-beta-D-galactopyranoside 
(Ci 2 FDG) (Debacq-Chainiaux et al. 2009). 
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Senescence induction by gamma-irradiation 

MSC were treated with a blood product irradiator (model IBL-437, 
Institute for transfusion medicine, Aachen) with 5, 10, 15, 20, or 25 
Gy as indicated in the text. Proliferation rate of MSC in early and 
late passage as well as after irradiation was determined in triplicate 
using the Alamar blue assay (Life Technologies), which was mea- 
sured with a Tecan infinite M200 plate reader. 

TERT overexpression 

The TERTsequence was excised from the pGRN145 plasmid (ATCC 
MBA- 141) and subcloned into the retroviral expression vector 
pMIG (pMSCV-IRES-GFP) (Addgene plasmid 9044, Addgene) via 
EcoRI. MSC of five donors were transduced at passage 4 with the 
TERT-IKES-GFP construct, or for control with the corresponding 
empty pMIG vector (IRES-GFP), and the culture expanded in par- 
allel until entering the senescent state. 

Conditional immortalization with TERT/SV40-TAg 

MSC were transfected with two lentiviral vectors (May et al. 2010): 
(1) SV40 large T antigen (SV40-TAg) under an unidirectional tet- 
racycline-dependent promoter, which also drives a reverse trans- 
activator through EMCV-IRES; (2) TERT under a bidirectional tet- 
racycline-dependent promoter, additionally driving expression of 
a reverse transactivator and an eGFP/neomycin phosphotransfer- 
ase fusion gene linked by an EMCV-IRES element. Initially, MSC 
were transduced on 2 consecutive d with SV40-TAg and after 
2 wk with TERT Upon neomycin selection, the resulting colonies 
(—20 independent clones) were pooled and culture expanded in 
DMEM (PAA; 4,5 g/L glucose) supplemented with glutamine (PAA), 
penicillin/streptomycin (PAA), 15% FCS (PAA), and doxycycline 
hyclate (Dox; final cone. 2 |xg/mL; Sigma, Taufkirchen). Proliferation 
rate with and without Dox was analyzed using the Alamar blue assay 
as described above. 



Telomere length analysis 

The average telomere length was estimated by monochrome mul- 
tiplex Q-PCR as described before (Cawthon 2009; Ziegler et al. 
2012). In brief, telomeres (T) and the beta-globin gene (S; single 
copy) were amplified on a MyIQ2 Two-color Realtime PCR De- 
tection-System (Biorad). Analysis was performed in triplicate in 
parallel to serial dilutions of standard DNA. The relative telomere 
to single-copy gene (T/S) ratio for experimental DNA samples was 
calculated by dividing the copy number of the telomere template T 
by the copy number of the beta-globin template S. Amplification 
efficiency for the telomere and the beta-globin was between 90% 
and 100% in each run. Specificity of the amplification was con- 
firmed by melt-curve analysis. 

Generation and characterization of iP-MSC 

MSC from three donors (passage 3) were infected on 2 consecu- 
tive d with equal amounts of pMXs-based retroviruses (Addgene) 
carrying the POU5F1, SOX2, KEF4, and MYC genes (Takahashi 
et al. 2007). On day 3, the medium was changed to high-glucose 
DMEM (Gibco), 20% FCS, 1% nonessential amino acids (NAA), lx 
penicillin/streptomycin, lx L-glutamine, 0.1 mM p-mercaptoethanol 
OME), and 50 ng/mL basic fibroblast growth factor (bFGF) (Peprotech) 
(Cai et al. 2010). On day 6 the cells were plated on gelatin-coated 
dishes containing irradiated CF1 murine embryonic fibroblasts 
(MEFs), and on day 7 the culture medium was supplemented with 
20 |xg/mL vitamin C (Sigma) and 1 mM valproic acid (Sigma). 



Between days 21 and 30 the ESC-like colonies were observed and 
mechanically isolated. Established iP-MSC were maintained on 
MEFs in DMEM/F12 medium supplemented with Glutamax, 20% 
knockout serum replacer, 1% NAA, 0.1 mM (3ME, and 50 ng/mL 
bFGF. Cells were passaged by manual dissection of cell clusters 
every 5-6 d. Human ESC lines HI, H9, and HES-2 (WiCell Research 
Institute) were used as controls and cultured under the same con- 
ditions as described above. Work with human ESC has been ap- 
proved by the Robert Koch Institute (permit no. Az 1 710-79-1-4-57). 
Immunohistochemical analysis of pluripotency markers was 
performed with primary antibodies against POU5F1 (Santa Cruz), 
SOX2 (Stemgent), NANOG (Stemgent), or TRA-1-80 (Santa Cruz) 
by alkaline phosphatase staining (Sigma) and by flow cytometric 
analysis of TRA-1-60 (BD Pharmingen) and SSEA4 (Santa Cruz). 
Furthermore, iP-MSC differentiation was analyzed using the em- 
bryoid body differentiation system with immunostaining after 
14-19 d for desmin (mesoderm), alpha-fetoprotein (endoderm), 
and nestin (ectoderm). Generation of iP-MSC clones and analysis of 
pluripotency has been described in further detail elsewhere (Shao 
et al. 2012). 

DNA methylation analysis 

Genomic DNA was isolated from 10 6 cells using the QIAGEN 
DNA Blood Midi-Kit. Bisulfite-converted DNA (EZ DNA Methyla- 
tion Kit; Zymo) was prepared and then analyzed using the 
InfiniumHumanMethylation450 BeadChip (Illumina) according to 
the manufacturer's instructions. DNAm is provided as beta-values 
for each of the 485.577 CpG sites, ranging from 0 (nonmethylated) 
to 1 (100% methylation). Hybridization and initial data analysis 
with the BeadStudio Methylation Module was performed at the 
DKFZ Gene Core Facility in Heidelberg. Distribution of samples 
across the BeadChips is presented in Supplemental Table 4. Data 
have been deposited at NCBIs Gene Expression Omnibus (GEO, 
http://www.ncbi.nlm.nih.gov/geo/; accession no: GSE37066). 

In addition, we analyzed all samples using our Epigenetic- 
Senescence-Signature, which has been described in detail before 
(Koch et al. 2012). This method is based on continuous DNAm 
changes at six specific CpG sites (associated with the genes GRM7, 
CASR, PRAMEF2, SEEP, CASP14, and KRTAP13-3). DNAm at these 
specific CpG sites was analyzed by pyrosequencing (Varionostic 
GmbH) and subsequently used for linear regression models to es- 
timate passage numbers, cPD, and days in culture. These calcula- 
tions can also be performed with the online calculator, which is 
accessible under http://www.molcell.rwth-aachen.de/dms/. 

Bioinformatic analysis of DNAm profiles 

For further analysis we have only considered CpG sites on au- 
tosomes. Initially, we considered various normalization regi- 
mens, but we then decided to use raw data beta-values due to the 
global changes of DNAm upon reprogramming with an overall 
increase of DNAm level. Histograms of DNA methylation levels 
and principal components analysis (PCA) were calculated with 
R. Unsupervised hierarchical clustering according to Euclidian 
Distance (with Ward's Minimum Variance Method) was performed 
using the MultiExperiment Viewer (MeV, TM4.7.4; http://www.tm4. 
org/mev/) or R. SA-DNAm changes were calculated by the limma 
package in R with adjusted P-values <0.05 and >20% change in 
DNAm level between early and late passages. Affiliation of CpG sites 
with gene regions or CGI was used as described in detail previously 
(Sandoval et al. 2011). Gene Ontology analysis of genes with DNAm 
changes was performed using the GoMiner software (http://discover. 
nci.nih.gov/gominer/) and enrichment of chromosomal location 
was analyzed with the gene set enrichment analysis (GSEA; http:// 
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www.broadinstitute.org/gsea/index.jsp). Furthermore, we compared 
DNAm profiles with published chromatin immunprecipitation 
(ChIP) results: data for H3K9me3 and H3K27me3 in adipose 
tissue-derived MSC were retrieved from GSE17053 (Delbarre et al. 
2010) and further processed by NimbleScan; EZH2 targets in bone 
marrow MSC were retrieved from the supplemental material of 
another manuscript (Wei et al. 2011). For comparison, DNAm and 
ChIP data were matched by GenBank Accession IDs and enrich- 
ment of SA-DNAm in these regions was determined by hyper- 
geometric distribution. 

Gene expression analysis 

For gene expression analysis, three additional MSC samples were 
isolated and expanded as described above. Reseeding was per- 
formed at a density of 13,333 cells/cm 2 in standard medium con- 
sisting of DMEM (PAA; 1 g/L glucose) supplemented with gluta- 
mine (PAA), penicillin/streptomycin (PAA), and 10% FSC (PAA). 
RNA was isolated at early passage (P3, P4, or P5) and late passage 
(P10) with the miRNeasy kit (Qiagen). RNA concentration and 
integrity were determined by an Agilent 2100 Bioanalyzer (Agilent 
Technologies, Inc.). Gene expression analysis was performed on 
GeneChip Human Gene 1.0 ST Arrays (Affymetrix) at the De- 
partment of Population Genetics, Genomics, and Bioinformatics 
at the University of Maastricht. Raw data were normalized by RMA 
using Affymetrix Power Tools, and matched to DNAm data by gene 
name. Data are accessible at GEO Series GSE37065. To provide 
further evidence that iP-MSC clones are fully reprogrammed we 
have analyzed gene expression profiles of iP-MSC clones as de- 
scribed above. These hybridizations were performed at the In- 
terdisciplinary Center for Clinical Research (IZKF) at the RWTH 
Aachen, and subsequently projected to PluriTest analysis as de- 
scribed before (Muller et al. 2011). All samples were clearly asso- 
ciated with pluripotent cells and data are accessible at GEO Series 
GSE38806. 

Data access 

The microarray data from this study have been submitted to the 
NCBI Gene Expression Omnibus (GEO) (http://www.ncbi.nlm. 
nih.gov/geo) under accession numbers GSE37065, GSE37066, and 
GSE38806. 
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